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(54) A region growing based noise reduction method for digital images 



(57) A method for removing noise from a pixel of a 
digital image channel includes the steps of: storing an 
original value of the pixel of interest and using the value 
of the pixel of interest and neighboring contiguous pix- 
els to determine a variable shape neighborhood region 



of cleaning pixels. The method further includes using 
the neighborhood region of cleaning pixels and the 
value of the pixel of interest to replace the original value 
of the pixel of interest with a noise cleaned pixel value. 
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Description 

[0001] The present invention relates to providing fil- 
tered digital images with reduced noise. 
[0002] Many noise reduction algorithms can be 
classified as non-linear spatial filters image processing 
algorithms. Most often these algorithms involve using 
the pixel values in a small local surrounding neighbor- 
hood to the pixel of interest combined with some form of 
non-linear weighting and/or statistical conditions 
applied to derive a noise free estimate. The small local 
surrounding neighborhood is usually centered on the 
pixel of interest. For this class of noise reduction algo- 
rithm the filter size is fixed, meaning that all image pixels 
are processed with the same size local surrounding 
neighborhood size. 

[0003] An example of a fixed size noise reduction 
algorithm is the Sigma Filter, described by Jon Son Lee 
in the journal article Digital Image Smoothing and the 
Sigma Filter, Computer Vision, Graphics, and Image 
Processing Vol 24, p. 255-269, 1983. This is a noise 
reduction fitter that uses a non-linear pixel averaging 
technique sampled from a rectangular window about 
the center pixel. Pixels in the local neighborhood are 
either included or excluded from the numerical average 
on the basis of the difference between the pixel and the 
center pixel. Mathematically, the Sigma Filter can be 
represented as 

qmn = £ ij a ijPij /r ij a ij 

and 

a ij = 1 if, Pij-Pmn' <=£ 
a ij = 0if i Pij -p mn t>e 

where py represents the pixels in the local surround 
about the center pixel p mn , q mn represents the noise 
cleaned pixel, and e represents a numerical constant 
usually set to two times the expected noise standard 
deviation. 

[0004] The Sigma Filter was designed for image 
processing applications for which the dominant noise 
source is Gaussian additive noise. Signal dependent 
noise sources can easily be incorporated by making the 
e parameter a function of the signal strength. However, 
for both signal independent and signal dependent noise 
cases the expected noise standard deviation must be 
known to obtain optimal results. The Sigma Filter per- 
forms well on highly structured areas due to the fact that 
most of the image pixels in the local neighborhood are 
excluded from the averaging process. This leaves high 
signal strength regions nearly unaltered. The filter also 
works well in large uniform areas devoid of image signal 
structure due to the fact that most of the local pixels are 
included in the averaging process. For these regions, 
the Sigma Filter behaves nearly as a low pass spatial fil- 



ter. 

[0005] Regions in images characterized by low 
amplitude signal modulation, or low signal strength, are 
not served well by the Sigma Filter. For these regions, 

5 most of the local pixel values are included in the averag- 
ing process thus resulting in a loss of signal modulation. 
Setting the threshold of the filter to a lower value does 
reduce the loss of signal, however, the noise is left 
mostly the same. 

70 [0006] Another example of a fixed size non-linear 
noise filter was reported by Arce and McLoughlin in the 
journal article Theoretical Analysis of the Max/Median 
Filter, IEEE Transactions Acoustical & Speech Signal; 
Processing, ASSP-35(1), p. 60-69, 1987 they named 

15 the Max/Median Filter. This filter separated the local 
4 surround region into four overlapping regions — hori- 
zontal, vertical, and two diagonal pixels with each 
region containing the center pixel. A pixel estimate was 
calculated for each region separately by applying a tak- 

20 ing the statistical median pixel value sampled from the 
regions' pixel values. Of these four pixel estimates, the 
maximum valued estimate was chosen as the noise 
cleaned pixel. Mathematically the Max/Median Filter 
can be represented as 

25 

qy = maximum of { Z 1f Z 2 , Z 3 , Z 4 } 

Zj = median of { Pij_ w , ...pjj,...,p ij+w } 
Z 2 = median of { p|_ w j, ... pj Jf ... t p i+w j } 
30 Z 3 = median of { Pj+wj.w, ... Pjj,...,Pj- W j +w } 
Z4 = median of { Pi_ Wi j_ W) ... Pi,j,...,Pj +Wl j+w 1 

Where qy represents the noise cleaned pixel, Z1, Z2, 
Z3, and Z4 represent the four pixel estimates, and py 

35 represents the local pixel values. The Max/Median Filter 
also reduces the noise present while preserving edges. 
For Gaussian additive noise, the statistical median 
value does not reduce the noise by as great a factor as 
numerical averaging. However, this filter does work well 

40 on non-Gaussian additive noise such as spurious noise. 
[0007] Noise is most visible and objectionable in 
images containing areas with little signal structure, for 
example blue sky regions with little or no clouds. The 
Sigma filter can produce a blotchy, or mottled, effect 

45 when applied image regions characterized by low signal 
content. This is largely due to the rectangular geometric 
sampling of local pixels strategy. The radial region sam- 
pling strategy employed by the Max/Median Filter pro- 
duces noise reduced images will less objectionable 

so artifacts in image regions characterized by low signal 
content. For images with high noise content, the arti- 
facts produced by radial region sampling strategy have 
a structured appearance. 

[0008] US-A-5,671 ,264 describes a variation of the 
55 Sigma Filter and Max/Median Filter. This algorithm bor- 
rows the technique of radial spatial sampling and multi- 
ple pixel estimates from the Max/Median Filter. 
However, the algorithm expands the number of radial 
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line segment to include configurations with more than 
four segments. The algorithm uses combinations of 
Sigma and Median filters to form the individual region 
pixel estimates. These pixel estimates derived from the 
N regions are then combined by numerical averaging or 5 
taking the statistical median value to form the noise 
cleaned pixel value. A key component of this algorithm 
is the randomization of one of the three essential region 
parameters: length, orientation, and number of regions. 
The randomization of the filter parameters is performed 10 
on a pixel to pixel basis thus changing the inherent char- 
acteristics with pixel location. It is claimed that the rand- 
omization feature reduces the induced structured 
artifacts produced by the radial region geometry sam- 
pling method. The algorithm described in US-A- 75 
5,671 ,264 can be categorized as a variable size non-lin- . 
ear noise filter. 

[0009] A. Lev, S.W. Zucker, and A. Rosenfeld, 
described two noise reduction algorithms in their journal 
article Iterative Enhancement of Noisy Images, IEEE 20 
Trans. Sysst. Man and Cybern. SCM-7, p. 435-441, 
1 977. Both these algorithms were based on edge sensi- 
tive local weighted averaging techniques. In the first 
algorithm, local pixel weights are assigned based on the 
presence or non-presence of edges produced with an 25 
edge detection algorithm. Pixels located on edges do 
not effect the local averaging process. Additionally, non- 
edge pixels separated from the center pixel by an edge 
pixel are also excluded (given a weight of zero). In this 
regard, only the region of non-edge pixels surrounding 30 
the center pixel is used to form the noise free pixel esti- 
mate. The second algorithm involved a generalization of 
the edge detection logic that formed the local pixel 
weights by combining the coefficients derived from four 
directional gradients filters, that is a horizontal, vertical 35 
and two diagonal gradient filters. The output of the 
directional gradient filters were transformed with an 
exponential function and combined in a multiplicative 
mode! to form the final 3 by 3 spatial filter. This filter is 
actually a filter of filtered responses. Both of the Lev et. 40 
al. filters reduce noise in images but must be applied 
iterative ly to have much effect due to small filter size (3 
x 3). Although these filters used a fixed size local sur- 
rounding neighborhood, they did incorporate the con- 
cept of confining the pixels of consideration to a local 45 
domain about the center pixel based on a detected fea- 
ture. However, the detection of that feature required the 
algorithm to have knowledge of pixels outside the local 
domain. 

[0010] Fixed local surround size noise reduction fil- so 
tens suffer from two main problems. Algorithms employ- 
ing small filter sizes take less computation time and 
preserve desirable low amplitude modulation signals 
but are also less affective at removing noise in unstruc- 
tured regions. Algorithms employing large filter sizes 55 
take more computation time and are more affective at 
removing noise in unstructured regions but also destroy 
desirable low amplitude modulation signals. Radial 



4 

region based noise reduction algorithms are affective at 
removing noise in unstructured regions but produce 
unwanted structured patterns in the noise cleaned 
images. Varying the size of the radial regions randomly 
can reduce the objectionability of the unwanted struc- 
tured patterns but does not impact the computation time 
required. 

[0011] It is an object of the present invention to pro- 
vide a noise reduction algorithm which combines the 
elements of variable filter size based on image pixel 
data with a radial regions. 

[0012] The present invention provides a noise 
reduction method that uses a variable shape and size 
filter by region growing a local neighborhood region of 
cleaning pixels, contiguous with the pixel of interest, 
based on the digital image pixel values. 
[0013] This object is achieved with a method for 
removing noise from a pixel of a digital image channel 
comprising the steps of: 

(a) storing an original value of the pixel of interest; 

(b) using the value of the pixel of interest and neigh- 
boring contiguous pixels to determine a variable 
shape neighborhood region of cleaning pixels; and 

(c) using the neighborhood region of cleaning pixels 
and the value of the pixel of interest to replace the 
original value of the pixel of interest with a noise 
cleaned pixel value. 

[0014] The present invention overcomes the limita- 
tion of fixed filter size algorithms by growing a region of 
contiguous pixels about the pixel of interest. Since the 
variable size and shape of the neighborhood region of 
cleaning pixels is based on the supplied pixel values, 
iow amplitude image detail is preserved while minimiz- 
ing the effect of induced structured artifacts. In addition, 
in highly structured regions of digital images, the 
present invention reduces the number of evaluated pix- 
els thus reducing the computation time. 

FIG. 1 is an example of the variable shape and size 

of the neighborhood region of cleaning pixels that 

can result from the application of the preferred 

embodiment of the present invention; 

FIG. 2 is an overall block diagram of the processing 

flow for practicing the present invention to produce 

digital images with reduced noise; 

FIG. 3 shows a flow chart of an algorithm which can 

accomplish all the processes of FIG. 2; 

FIG. 4 is a diagram of the radial path configuration 

used by the present invention; 

FIG. 5 is a diagram of the area fill propagation paths 

used by an alternative embodiment of the present 

invention; 

FIG. 6 is a diagram of the area fill path configuration 
used by an alternative embodiment of the present 
invention; and 

FIG. 7 is an example of the variable shape and size 
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of the neighborhood region of cleaning pixels that 
can result from the application of an alternative 
embodiment of the present invention. 

[0015] A digital image is comprised of one or more s 
digital image channels. Each digital image channel is 
comprised of a two-dimensional array of pixels. Each 
pixel value relates to the amount of light received by the 
imaging capture device corresponding to the physical 
region of pixel. For color imaging applications a digital 10 
image will often consist of red, green, and blue digital 
image channels. For monochrome applications, the dig- 
ital image will only contain one digital image channel. 
Motion imaging applications can be thought of as a 
sequence of digital images. Those skilled in the art will 75 
recognize that the present invention can be applied to, 
but is not limited to, a digital image channel for any of 
the above mentioned applications. In fact the present 
invention can be applied to any two dimensional array of 
noise corrupted data to obtain a noise cleaned output 20 
Although the present invention describes a digital image 
channel as a two dimensional array of pixels values 
arranged by rows and columns, those skilled in the art 
will recognize that the present invention can be applied 
to mosaic (non rectilinear) arrays with equal effect 25 
Those skilled in the art will also recognize that although 
the present invention describes replacing original pixel 
values with noise cleaned pixel values, it is also trivial to 
form a new digital image with the noise cleaned pixel 
values and retain the original pixel values in tact. 30 
[0016] FIG. 2 shows a functional block diagram of 
the present invention applied to a pixel of interest A 
pixel has two attributes that are used: a numerical value 
and a location within the two-dimensional array digital 
image channel which will be referred to as the pixel 35 
coordinates. Pixel coordinates include two indices 
which correspond to the row and column location. The 
original value of the pixel of interest, which will be 
referred to as the original pixel value, is denoted by X 
and is input to the statistics accumulator 30. The pixel 40 
coordinates of the pixel of interest, which will be referred 
to as the original pixel coordinates, are input to the pixel 
path selector 1 0. 

[0017] Each embodiment of the present invention 
contains one or more pixel processing paths each hav- 45 
ing a predetermined route of contiguous pixels. The 
pixel path selector 10 keeps track of which pixel 
processing path is currently being processed and sends 
an index indicating the current pixel processing path to 
the next pixel finder 20. The next pixel finder 20 deter- so 
mines the next sequential pixel to be evaluated. This is 
accomplished by retrieving the next sequential pixel of 
the current pixel processing path. The value of the next 
sequential pixel is fed to the statistics accumulator 30 
which uses it to determine if the next sequential pixel 55 
should be included in the neighborhood region of clean- 
ing pixels constituting a variable shape contiguous col- 
lection of pixels surrounding the pixel of interest. 



[0018] The processing of pixels by the statistics 
accumulator 30 is complete when the last pixel process- 
ing path is completed. Each pixel processing path is 
completed when one or more pixel path ending condi- 
tions are satisfied. Examples of pixel path ending condi- 
tions are: the location of the next sequential pixel lies 
outside the boundary of the two-dimensional array of 
pixels of the digital image channel, the total number of 
pixels processed by the statistics accumulator 30 has 
reached a maximum allowable number, the total 
number of pixels processed by the statistics accumula- 
tor 30 on the current pixel processing path has reached 
a maximum allowable number, and the noise cleaned 
pixel value output by the statistics accumulator 30 has 
not changed from its previous value by more than an 
allowable difference. When the statistics accumulator 
30 determines that a pixel path ending condition has 
been satisfied, a signal is sent to the pixel path selector 
1 0 to increment the index of the cu rrent pixel processing 
path. 

[001 9] When the processing of all the pixel process- 
ing paths is complete, the neighborhood region of 
cleaning pixels for the pixel of interest is defined. The 
neighborhood region of cleaning pixels is input to the 
noise cleaning calculator 40 which uses some or all of 
the neighborhood region of cleaning pixels to calculate 
a noise cleaned pixel value. The noise cleaned pixel 
value replaces the original value of the pixel of interest 
providing a pixel value with less inherent noise. This 
process is repeated for all the pixels in the digital image 
channel. Thus with the application of the present inven- 
tion the noise in the original digital image channel is 
removed. 

[0020] It is important to note that since each individ- 
ual pixel processing path is contiguous with the pixel of 
interest, so too is the neighborhood region of cleaning 
pixels. The present invention differs from all previous 
inventions in that the neighborhood region of cleaning 
pixels, which the noise cleaned pixel value is calculated 
from, is contiguous with the pixel of interest, variable in 
shape, with a shape which is completely determined by 
the pixel values contained in the original digital image 
channel. 

[0021] Referring to FIG. 3, a flow chart in block dia- 
gram form of an algorithm is shown which can practice 
the present invention as depicted in FIG. 2. The 
processing steps begin with block 1 00 which shows the 
step of finding a pixel of interest. The preferred embodi- 
ment of the present invention starts with the first pixel in 
computer memory at the first column and the first row. 
Block 1 1 0 shows the step of storing the original value of 
the pixel of interest. Block 120 shows the step of deter- 
mining a variable shaped neighborhood region of clean- 
ing pixels. This is performed by using the stored value of 
the pixel of interest and the values of pixel located in a 
local neighborhood about the pixel of interest. Block 130 
shows the processing step where a noise cleaned pixel 
value is calculated from the values of pixels included in 
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the variable shaped neighborhood region of cleaning 
pixels and the stored value of the pixel of interest with 
the application of noise reduction logic. Block 140 
shows the processing step of replacing the original 
value of the pixel of interest with the noise cleaned pixel 5 
value. The processing procedure then continues to 
block 100 where a new pixel of interest is found. The 
preferred embodiment of the present invention finds the 
next sequential pixel in memory as the next pixel of 
interest to be processed. The processing steps* shown w 
by blocks 100, 110, 120, 130, and 140 are repeated for 
each pixel of interest to be processed. The preferred 
embodiment of the present invention repeats these 
processing steps for all pixels in the digital image chan- 
nel and for all digital image channels in the digital 75 
image. 

[0022] Many possible pixel processing paths can be 
constructed. The preferred embodiment of the present 
invention employs a set of eight radially oriented pixel 
processing paths depicted in FIG. 4. In this regard, the 20 
eight pixel processing paths are disposed along radial 
directions in relation to the pixel of interest. This collec- 
tion of pixel processing paths will be called the radial 
path configuration. These eight pixel processing paths 
are labeled A, B, C, D, E, F, G, and H with sequential pix- 25 
els processed in each pixel processing path labeled 1, 
2, 3, and so on. Each pixel processing path starts at an 
adjacent pixel coordinate to the original pixel coordi- 
nate. As an example, path a shown in FIG. 4 starts by 
sending the pixel value corresponding to pixel location 30 
a1 to the statistics accumulator 30. If after processing 
this pixel value no pixel path ending condition is satis- 
fied, the next pixel processed will be located at a2. 
[0023] Since each pixel processing path must sat- 
isfy a pixel path ending condition and the pixel path end- 3S 
ing conditions can be based on the underlying pixel 
values encountered, the number of pixels processed in 
each pixel processing paths is allowed to vary. FIG. 5 
illustrates an example of a pixel region that can result 
using the radial path configuration. In this example the 40 
pixel path ending condition for each of the eight pixel 
processing paths a through h was satisfied at pixel loca- 
tions a6, b5, c4, d3, e3, f3, g4, and h4 respectively. The 
effective size and shape of the neighborhood region of 
cleaning pixels thus changes based on when the pixel 45 
path ending condition is satisfied for each individual 
pixel processing path. FIG. 1 illustrates the effective size 
and shape of the neighborhood region of cleaning pixels 
corresponding to the example shown in FIG. 4. The 
pixel marked X is the location of the pixel of interest. The so 
pixels marked Y represent the pixel locations a6, b5, c4, 
d3, e3, f3, g4, and h4 for which the pixel path ending 
condition was satisfied for each individual pixel process- 
ing path. 

[0024] An alternative embodiment of the present 55 
invention has a different set of pixel processing paths 
shown in FIG. 6. This set of four pixel processing paths 
will be called the area fill path configuration. Each pixel 



processing path starts out at a pixel location adjacent to 
the pixel of interest. The first pixel processing path des- 
ignated as path A, is defined by a series of vertical and 
horizontal steps. Starting with pixel a1 , the second pixel 
to be evaluated is a2 located above pixel a1. Since this 
pixel lies on a diagonal line emanating from the pixel of 
interest, the route of path A proceeds to the right to eval- 
uate pixel a3. Path A continues down until it reaches 
another pixel located on a diagonal line emanating from 
the pixel of interest. The other pixel processing paths B, 
C, and D follow the same logic but cover the pixel loca- 
tions above, to the left and below the pixel of interest 
respectively. FIG. 7 illustrates an example neighbor- 
hood region of cleaning pixels that can result with the 
application of the area fill path configuration. The pixel 
path ending condition was satisfied at pixels locations 
a8, b7, c21, and d5 for the individual pixel processing 
paths. 

[0025] The area fill path configuration does not 
leave any gaps, or missing pixel locations surrounding 
the pixel of interest while the radial path configuration 
does. Both of these embodiments constitute a region 
growing based noise filter. The effective size and shape 
of the noise filter depends on the choice of pixel 
processing paths but more importantly on the proc- 
essed pixel values. 

[0026] US-A-5,671 ,264 describes a noise filter 
which relies on a set of radially oriented regions consti- 
tuting a neighborhood region of cleaning pixels. A key 
aspect of this algorithm is based on the radomization of 
one of three attributes of the radially oriented neighbor- 
hood region of cleaning pixels: the number of regions, 
the length of the regions, and the orientation angle of 
the regions. By randomizing one or more of these 
attributes it is claimed in US-A-5,671 ,264 that the result- 
ing noise cleaned image is improved. The method 
described in US-A-5,671 ,264 does constitute a variable 
size and shape noise filter. However, the variable size 
and shape of the noise filter is based not of the digital 
image pixel values but by a random number generator. 
[0027] The algorithms described by A. Lev, S.W. 
Zucker, and A. Rosenfeid, in their journal article Iterative 
Enhancement of Noisy Images, IEEE Trans. Sysst. Man 
and Cybern. SCM-7, p. 435-441, 1977 can also be con- 
sidered to be a variable size noise filter. In these algo- 
rithms some of the pixels are excluded from contributing 
to the calculation of the noise cleaned pixel value based 
on presence of detected edge information. However, in 
Lev et af's algorithm all of the pixel values within a fixed 
size and shape filter region must be evaluated before 
the excluded pixels can be identified. In this regard, the 
Lev et. al algorithm can be considered a fixed size noise 
filter with variable pixel value weights. 
[0028] The preferred embodiment employs a pixel 
path ending condition, used by the statistics accumula- 
tor 30, and a method for calculating the noise cleaned 
pixel value based, used by the noise cleaning calculator 
40, based on the Sigma Filter, described by Jon Son 
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Lee in the journal article Digital image Smoothing and 
the Sigma Fitter, Computer Vision, Graphics, and 
Image Processing Vol 24, p. 255-269, 1983. Pixel val- 
ues sent to the statistics accumulator 30 shown in FIG. 
2 are compared with the value of the pixel of interest. If 5 
the difference between the value of the next sequential 
pixel and the value of the pixel of interest is less than a 
threshold e, the next sequential pixel is included in the 
neighborhood region of cleaning pixels and by the noise 
cleaning calculator to calculate a numerical average io 
which forms the noise cleaned pixel value. The numeri- 
cal constant e is usually set to two times the expected 
noise standard deviation. Mathematically, the numerical 
average can be expressed as 

15 

c lmn = E ij a jk Pjk 1 Z jk a jk 

and 

a jk = 1 if lp jk -p mn l <=e 20 

a jk = 0if 'Pjk-Pmn'^ 

where pj k represents the k th sequential pixel in the j m 
pixel processing path, p mn represents the value of the 25 
pixel of interest located at row m and column n, and q mn 
represents the noise cleaned pixel value. Thus the 
noise cleaned pixel value is based on the difference 
between the value of the pixel of interest and the values 
of pixels and is represented by the statistical mean of 30 
the values of the pixels contained in the neighborhood 
of cleaning pixels. The pixel path ending condition for 
each pixel processing path is based on the number N of 
excluded pixels (pixels for which aj k equals zero) in the 
individual pixel processing path expressed as 35 

Af = k-£ |c a jk . 

The variable k represents the number of processed next 
sequential pixel values encountered in the current pixel 40 
processing path. When the number of excluded pixels N 
is equal to a predetermined threshold T, the pixel path 
ending condition for the pixel processing path is satis- 
fied. The preferred embodiment of the present invention 
uses a value of two for the predetermined threshold T. 45 
This pixel processing path ending condition thus deter- 
mines the length of the predetermined pixel processing 
path, including paths disposed along radial directions, 
based on the difference of the value of the pixel of inter- 
est and the values of the pixels contained in the pixel so 
processing path. 

[0029] It should be noted that while the preferred 
embodiment of present invention uses the same numer- 
ical value for e to determine the noise cleaned pixel 
value q mn and the pixel path ending condition, this does ss 
not have to be the case. The preferred embodiment 
uses the same value for e for simplicity and computa- 
tional efficiency. If the source of noise is signal depend- 



ent, that is depends on the value of the pixel of interest 
Pmn. the value of e should track with p mn . An alternative 
embodiment of the present invention uses a look-up- 
table of e values. As each pixel of interest is processed, 
the value of p mn is used as an index into the look-up- 
table of e values to set the current value of e. 
[0030] The preferred embodiment also uses the 
boundary of the digital image channel as a pixel path 
ending condition for the pixel processing paths. This 
must be used since there is no pixel data beyond the 
boundary. Another pixel path ending condition used by 
the preferred embodiment is a maximum number of pix- 
els for each pixel processing path. Typically a maximum 
of 15 pixels for each pixel processing path is used. 
[0031] The pixel inclusion/exclusion logic of Lee's 
Sigma Filter can be represented mathematically as 

qmn^ijayPij^ij 3 ! 

and 

a ij = 1 if 'Pij-Pmn' <=* 
ajj = 0 if l Pij - Prnn l>e 

where p^ represents the pixels in the local surround 
about the original pixel p mn , q mn represents the noise 
cleaned pixel, and e represents a numerical constant 
usually set to two times the expected noise standard 
deviation. The Sigma Filter as described by Lee was 
applied to a fixed size rectangular pixel region centered 
about the original pixel. 

[0032] An alternative embodiment of the present 
invention uses a median filter to calculate the noise 
cleaned pixel value produced by the noise cleaning cal- 
culator 40. In this embodiment, the same pixel path end- 
ing conditions are used as in the preferred embodiment. 
Therefore the neighborhood region of cleaning pixels is 
determined by the Sigma filter mathematics. The noise 
cleaned pixel value is given by the statistical median of 
the pixel values contained in the neighborhood region of 
cleaning pixels. 

[0033] The present invention can be provided in a 
computer program which is stored on a computer read- 
able storage medium. Such a medium can comprise for 
example; a magnetic disk (such as a floppy disk), mag- 
netic tape, code bars, solid state electronic storage 
devices (such as random access memories or read only 
memories), or any other physical device or medium 
which can be employed to store a computer program. 

Claims 

1 . A method for removing noise from a pixel of a digital 
image channel comprising the steps of: 

(a) storing an original value of the pixel of inter- 
est; 
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2. 



3. 



(b) using the value of the pixel of interest and 
neighboring contiguous pixels to determine a 
variable shape neighborhood region of clean- 
ing pixels; and 

(c) using the neighborhood region of cleaning 
pixels and the value of the pixel of interest to 
replace the original value of the pixel of interest 
with a noise cleaned pixel value. 

A method for removing noise on a pixel by pixel 
basis from pixels of a digital image channel com- 
prising the steps of: 

(a) storing an original value of the pixel of inter- 
est; 

(b) using the value of the pixel of interest and 
neighboring contiguous pixels to determine a 
variable shape neighborhood region of clean- 
ing pixels; 

(c) using the neighborhood region of cleaning 
pixels and the value of the pixel of interest to 
change the original value of the pixel of interest 
so that it has been noise cleaned; and 

(d) repeating steps (a)-(c) for other pixels of 
interest. 



tained in the neighborhood of cleaning pixels. 

9. A computer program product, comprising; having a 
computer readable storage medium having a com- 
5 puter program stored thereon for removing noise 

from a pixel of a digital image channel comprising 
the steps of: 

(a) storing an original value of the pixel of inter- 
10 est; 

(b) using the value of the pixel of interest and 
neighboring contiguous pixels to determine a 
variable shape neighborhood region of clean- 
ing pixels; and 

15 (c) using the neighborhood region of cleaning 

pixels and the value of the pixel of interest to 
replace the original value of the pixel of interest 
with a noise cleaned pixel value. 

20 10. A computer program product, comprising; having a 
computer readable storage medium having a com- 
puter program stored thereon for removing noise on 
a pixel by pixel basts from pixels of a digital image 
channel comprising the steps of: 



25 



The method of ctaim 1 wherein step (b) the neigh- 
borhood region of cleaning pixels is chosen from 
one or more predetermined pixel processing paths 
contiguous with the pixel of interest. 30 

The method of claim 3 wherein step (b) the prede- 
termined pixel processing paths contiguous with 
the pixel of interest are disposed along radial direc- 
tions in relation to the pixel of interest. 35 

The method of claim 3 wherein step (b) further 
comprises determining the length of the pixel 
processing paths based on the difference of the 
value of the pixel of interest and the values of the 40 
pixels contained in the pixel processing paths. 



(a) storing an original value of the pixel of inter- 
est; 

(b) using the value of the pixel of interest and 
neighboring contiguous pixels to determine a 
variable shape neighborhood region of clean- 
ing pixels; 

(c) using the neighborhood region of cleaning 
pixels and the value of the pixel of interest to 
change the original value of the pixel of interest 
so that it has been noise cleaned; and 

(d) repeating steps (a)-(c) for other pixels of 
interest. 



6. The method of claim 4 wherein step (b) further 
comprises determining the length of the pixel 
processing paths based on the difference of the 45 
value of the pixel of interest and the values of the 
pixels contained in the pixel processing paths. 

7. The method of claim 1 wherein step (c) further com- 
prises calculating a noise cleaned pixel value so 
based on the difference between the value of the 
pixel of interest and the values of the pixels con- 
tained in the neighborhood of cleaning pixels. 

8. The method of claim 2 wherein step (c) further com- 55 
prises calculating a noise cleaned pixel value 
based on the difference between the value of the 
pixel of interest and the values of the pixels con- 
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